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ABSTRACT 

This paper presents further results from our spectroscopic study of the globular cluster 
(GC) system of the group elliptical NGC 3923. From observations made with the 
GMOS instrument on the Gemini South telescope, an additional 50 GC and Ultra 
Compact Dwarf (UCD) candidates have been spectroscopically confirmed as members 
of the NGC 3923 system. When the recessional velocities of these GCs are combined 
with the 29 GC velocities reported previously, a total sample of 79 GC /UCD velocities 
is produced. This sample extends to over 6 arcmin (>6 R e ^30 kpc) from the centre 
of NGC 3923, and is used to study the dynamics of the GC system and the dark 
matter content of NGC 3923. It is found that the GC system of NGC 3923 displays no 
appreciable rotation, and that the projected velocity dispersion is constant with radius 
within the uncertainties. The velocity dispersion profiles of the integrated light and 
GC system of NGC 3923 are indistinguishable over the region in which they overlap. 
We find some evidence that the diffuse light and GCs of NGC 3923 have radially 
biased orbits within ^130". The application of axisymmetric orbit-based models to 
the GC and integrated light velocity dispersion profiles demonstrates that a significant 
increase in the mass-to- light ratio (from M/Ly = 8 to 26) at large galactocentric radii 
is required to explain this observation. We therefore confirm the presence of a dark 
matter halo in NGC 3923. We find that dark matter comprises 17.5^4 5% of the mass 
within 1 Re, 41,2±^;§% within 2 R c , and 75.6±^;|% within the radius of our last 
kinematic tracer at 6.9 R . The total dynamical mass within this radius is found to 
be 1. 

5^0.25 xlO 12 M Q . In common with other studies of large ellipticals, we find that 
our derived dynamical mass profile is consistently higher than that derived by X-ray 
observations, by a factor of around 2. 

Key words: galaxies: general - galaxies: individual: NGC 3923 - globular clusters: 
general 
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The mass distributions and dark matter (DM) content of 
galaxies are of significant importance to the study of the 
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formation and evolution of galaxies. The current cold dark 
matter (CDM) paradigm posits that dark matter is the dom- 
inant matter component found in the Universe, and hence 
drives the formation of structure throughout the history of 
the Universe. However making measurements of the actual 
distribution of DM with which to confront theory can prove 
to be difficult. In particular, measuring the DM distribu- 
tion on the scales of individual galaxies can prove to be 
extremely trying for any galaxy type other than the most 
massive early-type galaxies (where significant X-ray emis- 
sion is present) or later-type spirals (where gas emission 
can be traced to larg e radii to trace flat rotation curves, 
e.g. IPersic et al.lll996l ). Intermediate to low mass early-type 
galaxies are of notable concern, as these galaxies usually 
lack significant X-ray or ionised gas emission, and gener- 
ally lack readily observable tracers of the stellar kinematics 
at large radii (such as absorption line kinematics), due to 
the inherent faintness of the stellar populations at the radii 
where DM is expected to dominate (beyond 2R e ). For these 
galaxy types then, the use of dynamical probes of the galaxy 
potential such as GCs and planetary nebulae (PNe) are of 
particular importance. 

GCs are well suited to this type of study for several rea- 
sons: (i) their ubiqui t y arou nd all relatively massive galaxies 
(e.g. see IPeng et ail (|2006l ) for Virgo Cluster galaxies); (ii) 
their extended distribution, that ensures significant numbers 
can be found beyond the radii at which the DM component 
becomes significant; (iii) and the fact that they are suffi- 
ciently luminous to be spectroscopically studied to large dis- 
tances (~30Mpc). In addition to providing determinations of 
the relative amounts of dark and baryonic matter, the study 
of the kinematics of large samples of GCs around early-type 
galaxies also provides other useful clues to their formation 
histories. In some galaxies for example, the GC system (or 
a subpopulation) is found to rotate with the galaxy (see be- 
low) , perhaps indicative of a common formation history for 
the stellar component and the GC system. 

The examination of the kinematics of GC systems is 
therefore an extremely promising area of research, with 
much to be learned about the processes that lead to the 
formation of currently observed galaxies. In recent years, 
a growing number of studies have made use of 8m class 
telescopes to build large samples of GC velocities around 
several early-type galaxies. In general, these have been of 
more massive cluster galaxies, known to harbour large GC 
populations. The project of which this study is a part in- 
tends to widen the host galaxy mass and environment r ange 
so far studie d. To date, NGC 33 79 (|Pierce et al.ll2006h and 
NGC 4649 (|Bridges et al.l 120061 ) have been studied. Both 
galaxies displayed GC velocity dispersion profiles that were 
constant with radius, a finding indicative of the presence of 
DM. 

To these studies c an be added investigations of the GC 
syste ms of NGC 139 9 (iRichtler et alll^jolfr ISchuberth et all 
|2010D, NGC 1407 ilRomanowskv et al 1 12009|) , NGC 447 2 



( Sharpies et~aH Il998l; jZepf et all |2000|; ICote et all l2003h . 
NGC 4486 (ICote et al ll200ll;IStrader et al.ll201ll), NGC 4636 



(ISchuberth et al J2006l ;1 Chakrabartv fc Ravchaudhury|2008l : 



IPark et alJl201Ct ). and NGC 5128 (jWoodlev et al.ll2007l) all 

of which have found similarly constant or slightly increasing 
GC velocity dispersion profiles, providing strong evidence for 
the existence of significant amounts of DM in these galaxies. 



Examination of the GC systems of these galaxies for 
signatures of rotation like that displayed by the Milky Way 
metal-rich GC population has been less conclusive. Sig- 
nificant rot ation has been observed in some galaxies (e.g. 



NGC 3115; 
NGC 5128; 
ers (e.g. NGC 1399 



Kuntschner et alj|2002l ; I Arnold et al] l201ll and 
Woodlev et al.l 120071) but little or none in oth- 



Richtler et ej"]|200Sl fc ISchuberth et al.1 



20ld. NGC 1407 ; iRomanowskv et al.l I2OO9I. NGC 3379 



Pierce et alj|200d . NGC 4472; IZepf et al.ll2000l ; ICote et all 
20031. NGC 4594; iBridges et al.l 120071 . and NGC 4649; 



Bridges et alll2006l ). In the cas e of NGC 4486 e arly studies 



indicated signifi cant rotation ( Coteet al.ll200ll ). but more 



recent analyses (|Strader et al.l 



201 ll ) indicate generally in- 



significant amounts of rotation. When GCs are split into 
blue and red subpopulations the picture becomes even more 
complicated. Generally it can be stated that the blue sub- 
population displays a higher velocity dispersion than the 
red population, but beyond this the picture is confused with 
considerable variation from galaxy to galaxy. 

The study of the orbits of the GCs has generally been 
stymied by th e need for extremely large s amples of GC veloc- 
ities (see e.g. iMerritt fc Tremblavlll99"3l . who estimate that 
around a thousand GC velocities would be required to solve 
for both the mass distribution and orbital anisotropy simul- 
taneously in the case of NGC 4486). Where fewer velocities 
are available, other input such as an X-ray profile can be 
used to constrain the potential, allowing the observed dis- 
persion profile to be used to determine the GC orbits. In 
NGC 1399 and NGC 4472 the GCs displ ay orbital character- 
istics close to isotropy, but in NGC 1407 (jRomanowskv et al.l 
120091 ) and NGC 4649 l|Bridges et al.ll200d ) there is some ev- 
idence of the GCs having more tangential orbits at certain 
radii. 

This paper exten ds the study of th e GC system of 
NGC 3923 described in lNorris et~all (|200Sh to include an in- 
vestigation of the kinematics of the NGC 3923 GC system. 
NGC 3923 is a large (R c of ~53.4 arcsec, 4. 6 kpc), nearby 
(D ~ 21 Mpc) elliptical galaxy (e ~ 0.3 - 0. 4. ISikkema et al. 
20071) . with a prominent shell structure (|Malin fc Carter! 
19831 ). and is the brightest galaxy of an average sized group. 
To date no studies have examined the GC kinematics of 
this galaxy, but several have pr obed the kinematic s of the 
integrated light of NGC 3923. In lNorris et~al] (|200Sl) diffuse 
light kinematics extracted from the same MOS slitlets as 
the GC spectra are used to confirm that the main body of 
NGC 3923 rotates along the major axis with a small ampli- 
tude of ^30 kms" 1 . This res ult is consistent with the work of 
iKoprolin fc Zeilingerl (|2000|). bu t in m arginal disagreement 



with the work of Carter et al.l (|l998l ) that found no rota- 
tion around the major axis but slight rotation of amplitude 
~20 kms -1 around the minor axis. The only previous inves- 
tigation to hav e studied the DM cont ent of NGC 3923 is the 
X-ray study of iFukazawa et all (|2006l ) , which found that the 
M/L in the B-band increased from 3.5 in the inner regions 
to around 15 at 18.4 kpc, a trend strongly supporting the 
presence of a DM halo associated with this galaxy. 



2 OBSERVATIONS AND DATA REDUCTION 

In this paper we continue t he study of the GC system of 
NGC 3923 first presented in lNorris et~aH (|2008h . To the 29 
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NGC 3923 GC velocities previously determined from Multi- 
Object Spectroscopy (MOS) we add an additional 50 GC 
velocities measured using a combination of MOS and Nod- 
and-Shuffle (N&S) observations. 

As described in lNorris et al.l (|2008l ), pre-imaging for GC 
candidate selection was undertaken for 3 fields (central, SW 
and NE) on 2004 January 19, comprising 4x200 seconds in 
Sloan g, and 4x100 seconds in r and i. Reduction of the 
imaging data involved bias subtraction, flat fielding, fringe 
removal (for i images), combining dithered sub-exposures 
to remove chip gaps, and finally removal of a median filter 
image (to remove the galaxy diffuse light). More details of 
the imaging r educt ion an d object select i on ca n be found in 
iForbes et all (|2004l ) and iBridges et all i|2006l ), with a full 
analysis of the imagi ng data from all ga laxies studied in this 
project presented in lFaifer et al.l (|201ll ). 

A full description of the reduction and analysis of the 
de ep MOS obs e rvatio ns of the central pointing is provided 
in lNorris et alj (|2008l ). We present here data from an addi- 
tional shallower (~3hrs vs. 8hrs) MOS observation centred 
on NGC 3923 and three further N&S observations. The shal- 
lower MOS observation was primarily focussed on studying 
the internal kinematics and stellar population s of a sample 
of UCDs discovered by iNorris fe Kannappanl (120111 ) . This 
MOS mask contained 33 slitlets, with those on UCD candi- 
dates being 0.5" wide and those on the fainter GC candidates 
1.5" wide. The B1200 grating was employed providing reso- 
lution of ~ 1.2 and 3.6 A and coverage of 1490 A generally 
centred around 5000A (the coverage varies slightly based on 
the position of the slit on the mask). These newer MOS ob- 
servations we re reduced us i ng the same tasks and procedures 
as outlined in lNorris et al.l (|2008l ). and interested readers are 
referred there for further details. Hereafter we shall focus on 
describing the reduction of the N&S spectroscopy. 

GMOS N&S masks were produced for each of the 3 
fields of pre-imaging. The central field containing the largest 
number of targets was therefore studied three times, twice 
with MOS spectroscopy and once with N&S. The N&S 
masks contained a total of 44 (central field), 46 (SW field) 
and 48 (NE field) slitlets, each 1 arcsecond wide by at least 
2.25 arcseconds long. The majority of the slitlets per mask 
(30-35) were placed on GC candidates, with the remaining 
slitlets placed on other objects of interest to fill up the mask. 
Each N&S mask was exposed using the B600_G5303 grism 
for 4x1800 seconds at a central wavelength of 505 nm and 
4x1800 seconds at a central wavelength of 510 nm, yielding 
4 hours of on-source integration per mask. The spectra were 
exposed through the g_G0325 and GG455_G0329 filters to 
restrict the wavelength coverage to a well-defined bandpass 
of ~4600-5550A. This restriction in wavelength coverage al- 
lowed more freedom in the positioning of slitlets and hence 
in choice of targets, as spectra could be aligned in the spatial 
dimension without overlapping. The wavelength coverage of 
each spectrum was still sufficient to cover most of the strong 
absorption lines found in optical spectra and includes the 
most important lines for such studies, H/3, Mg6, Fe5270 and 
Fe5335. Bias frames, flat fields and copper-argon (CuAr) arc 
spectra were observed throughout the observations as part 
of the standard Gemini baseline calibrations. 

The reduction of the N&S spectra was accomplished us- 
ing the standard Gemini/GMOS packages in 1RAF, a brief 
description of which shall be given here to help illustrate the 



benefits and problems associated with using this technique 
for GC studies. In N&S spectroscopy the target is placed 
at one end of the slit, it is exposed as normal for a period 
of time, and the telescope is then nodded to a sky position 
whilst shuffling the charge on the CCD between science and 
un-illuminated storage regions. The image that results con- 
tains two spectra - one of the object and one of the sky, both 
of which importantly were obtained on exactly the same pix- 
els after traversing identical optical paths. Therefore, when 
subtracting the sky spectrum from the object spectrum, ef- 
fects such as flat-fielding errors, fringing, and variation in 
the sky cancel out more accurately than in standard sky 
subtraction procedures. In the case of the data here, a small 
offset was applied in the nodding procedure so that the ob- 
ject appears in both the science and the sky spectra, only at 
different ends of the slit. To produce these spectra, the stan- 
dard reductions for Gemini/GMOS IRAF N&S tasks were 
used to provide bias subtraction, flat fielding, extraction of 
each spectra pair into separate extensions, and wavelength 
calibration. At this point it is possible to shift one spec- 
trum up to overlap with the other, so that the sky region 
of one spectrum is over the target spectrum of the other. 
Subtracting one spectrum from the other then produces a 
sky subtracted target spectrum pair, where one spectrum is 
negative and the other positive. It is then a simple matter 
to trace and extract each spectrum as normal and combine 
each exposure together (remembering to flip the negative 
spectra to positive ones) to produce the final ID spectrum. 

In theory, the use of the N&S method should have one 
major advantage over standard MOS observations: the abil- 
ity to pack more objects onto each mask. This is because 
the improved sky-subtraction accuracy relative to MOS ob- 
servations means that shorter slitlets are required to deter- 
mine equivalently accurate sky spectra. In the specific case of 
the data presented here, this increase in efficiency was not 
achieved due to significant issues with the GMOS instru- 
ment in N&S mode. In particular the GMOS CCDs exhibit 
small regions that do not allow the charge to move freely 
through them during the "shuffling" part of the N&S pro- 
cedure. These "charge traps" cause trails of charge deficit. 
The net effect of these charge traps is that they lead to 
artifacts in the final extracted ID spectra that are indistin- 
guishable from absorption features. Several methods have 
been developed to d eal with these traps (see for example 
I Abraham et alj|2004l ). However, no method is straightfor- 
wardly applicable to the data available here, as they gener- 
ally require additional steps during the observations. Hence, 
these charge traps make some regions of the GMOS CCD 
unusable in the present case, leading to the unfortunate loss 
of up to 25% (~10 objects) of the spectra per mask. 



2.1 GC Velocity Determination 

As in INorris et al.l (|2008l ). the recessional velocities of the 
GC candidates were measured using the FXCOR task in 
the RV package of IRAF. This task carries out a fourier 
cross correlation between the object spectra and a spec- 
tral library. In this case no radial velocity standard stars 
were observed as part of this project, therefore theoretical 
te mplate spectra f rom the simple stellar population library 
of IVazdekisI (1 19991) were used as templates. As descri bed in 



INorris et al.l ( 20081 ) a further 6 stellar spectra from the lJonesl 



© 2007 RAS, MNRAS OOO-fiHTol 



4 Norris et al. 



2500 a 




500 1000 1500 2000 2500 
Deep MOS Velocity [kms 1 ] 

Figure 1. Comparison of the velocities of 24 objects observed in 
the deep MOS mask and then re-observed on at least one of the 
shallow MOS or Nod-and-Shufnc fields. The dashed blue line is 
the best fit least squares linear relation which has the form Y = 
Xx 1.003 (± 0.017) - 9.79 (±26.62). The RMS scatter about the 
one-to-one relation is 62.3 kms -1 . 



(|l997h library were added to ensure adequate coverage of the 
low-metallicity region. Hence the template spectra spanned 
a range in metallicity from [Fe/H] < -1.5 to +0.2 and in 
age from 1-18 Gyr. The quoted velocity for each object is 
the 3-ct clipped mean of the velocities derived from the FX- 
COR fits to each template; the errors are estimated from 
the mean of the errors measured by FXCOR for those ve- 
locities not clipped. Only GC candidates where at least half 
of the templates returned Tonry and Davis R value greater 
than 5 and velocities within 3-ct of the mean are consid- 
ered further. This ensures that only reliable velocities are 
included in our analysis, and is a necessary step since only 
a few erroneous velocities can lead to significant changes in 
the inferred propert ies such as velocity dispersion (see e.g. 
IStrader et aLll201lf ). 

Objects with velocities in the range 1100-2500kms _1 
(~V g ai ± 700 kms -1 , or ~3-<r of the GC velocity disper- 
sion) are assumed to be associated wit h NGC 3923. This i s 
marginally different to that quoted in iNorris et al.l <|200Sh . 
where 1200-2400 kms -1 was used. In practice, however, this 
makes no difference t o the selection of objects classified as 
GCs in INorris et all (|2008l ). as no additional objects from 
the deep MOS observations would fall into this wider veloc- 
ity range. The wider velocity range does however allow one 
additional N&S object to be included in the present sample. 

A total of 24 objects were examined with both the deep 
MOS and shallow MOS or N&S allowing an examination of 
the robustness of the measured velocities and their errors. 
Figure [1] displays the velocity measured using the shallow 
MOS or N&S observations, versus the velocity measured us- 
ing the deep MOS observations for 24 objects that were ob- 
served twice. As can be seen, the agreement between the two 
measurements is exceptionally good, with the best fit linear 
relation between the two sets of observations being consis- 
tent with no offset between the two methods. This good 



agreement allows the simple weighted averaging of veloci- 
ties measured using both techniques in the case of overlaps. 

In the few cases where N&S targets were observed on 
two separate N&S masks but no velocity was recovered, the 
spectra were coadded and then velocity measurement reat- 
tempted. However, in no case were any additional velocities 
recovered. In total, after combination of repeat observations, 
79 secure velocities were determined. The spatial distribu- 
tion of these 79 GCs and UCDs, separated into blue and red 
subpopulations as described in Section 13.11 can be seen in 
Figure [2] Table [1] provides the positions, photometry (from 
IFaifer et aill201ll ) and velocities for the 79 objects confirmed 
as NGC 3923 GCs. 



3 RESULTS 
3.1 GC Colours 

In this subsection we examine the photometry of the 
NGC 3923 GCs in order to be able to robustly separate the 
GCs into blue and red subpopulations. Figure [3] displays the 
colour-magnitude diagram for all candidate GCs/UCDs de- 
tected in th e GM OS preimaging (grey dots), as found by 
IFaifer et al.l (|201ll ). Additionally, the 79 spectroscopically 
confirmed GCs/UCDs are overplotted as large filled orange 
circles/stars. In the upper panel of Figure [3] the respec- 
tive histograms of g-i colour for each of the two samples 
(all potential GCs/UCDs and spectroscopically confirmed 
GCs/UCDs) are also displayed. 

It can be seen from Figure [3] that a significant fraction 
of the brightest GC candidates were not spectroscopically 
observed. This is largely due to problems placing pairs of 
slitlets because of potential spectral overlaps. Other poten- 
tial GC candidates could not be observed or were subse- 
quently lost because the GCs were located too close to the 
galaxy centre, or because the spectra were unusable due to 
the previously described problems with the N&S observa- 
tions. 

The histogram of the full GC candidate sample displays 
a dip at around g -i ~ 0.87 characteristic of the existence of 
a bimodal colour distributio n. In their c o mpreh ensive anal- 
ysis of the GC photometry IFaifer et al.l l|201ll ) find unam- 
biguous bimodality in the colours of the GCs of NGC 3923 
(with a dividing colour of go-io = 0.89), in line with other 
studi es of NGC 3923 using bot h ground based dZepf et al. 
19951) , and HST photometry dSikkema et all 12006] ; IChc 



20081; INorris fc Kannappanl 1201 ll ). Importantly however, 



Faifer et al.l (|201ll ) find that bimodality is strongest for the 
fainter GCs and disappears for the brightest GCs. This 
result can be interpreted as evidence f or a blue til t (e.; 



Strader et all 120061; iMieske et all 120061 ; lHarris et all 120061 ; 
Peng et all 120091 ; iForbes et all 120101 1 in the colours of the 



blue GC population of NGC 3923, where the disappearing 
bimodality is simply due to the merging of the blue GCs 
into the red locus at high luminosities. 

The presence or absence of a blue tilt in the NGC 3923 
GC system is important. As it is only possible to measure ve- 
locities for the most luminous GCs (see Figure [3J, significant 
numbers of "red" GCs may actually have formed as part of 
the blue subpopulation. Therefore, studies of the kinemat- 
ics of GC subpopulations may suffer significant bias in their 
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Figure 2. The positions of kinematically confirmed GCs/UCDs associated with NGC 3923. Blue circles and red squares indicate blue 
and red GC sub-popu l ation members, respectively. The three stars indicate the positions of the three confirmed UCDs discovered by 
iNorris & Kannappanl j201ll ). The subpopulations were split at gp-in = 0.89. The phy sical scale is indicated based on our adopted 
distance for NGC 3923 of 17.9 Mpc (chosen to match that used in lFukazawa et aHl2006l ). 



derived properties, due to contamination of the red popula- 
tion by members of the blue population. Unfortunately, with 
the relatively limited dataset on hand, it is not possible to 
carry out a more sophisticated separation of blue and red 
population GCs (perhaps by allowing a sloping separation 
in colour-magnitude at higher luminosity) . Consequently, in 
all analyses to follow, where necessary the GCs/UCDs are 
separated using g-i = 0.89. This separation divides our total 
sample of 79 objects into subsamples of 34 blue and 45 red 
objects. 



3.2 GC Spatial Distributions 

The spatial distributions of GC subpopulations provide use- 
ful clues to their formation, as well as being essential inputs 
into the dynamical models to be described in Section [3. 61 

Figure [4] displays the GC surface density profiles of the 
total, blue, and red samples for the full photometrically se- 



lected sample from iFaifer et all (j201lh . The profiles were 
constructed by normalising GC number counts from HST 
ACS imaging of the inner regions (because it had higher 
completeness than the GMOS imaging in the inner region) 
to the wider field GMOS imaging in an overlap region of 
radius 60 to 90 arcsecs. The profiles were determined in cir- 
cular annuli, as no significant evidence of ellipticity in the 
GC system was found. They were also corrected for areal 
incompleteness, to take into account the arrangement of the 
GMOS imaging pointings. In each panel of Figure [4] the 
profiles have been normalised to match the NGC 3923 V- 
band surface brightness profile at 60 arcsec, allowing direct 
comparison between the GC and NGC 3923 diffuse light 
distributions. 

It can be seen that the red GCs are more centrally con- 
centrated and have a steeper number density profile than 
that of the blue GCs, a result common to many studies of 
the GC systems of early-type galaxies (see e.g. I Lee et aT] 



© 2007 RAS, MNRAS 00Q.[TJfT5l 



6 Norris et al. 




19- 



it 




0.4 0.6 0.8 1.0 1.2 1.4 
(g'-i')o [mag] 

Figure 3. Lower Panel: Colour-magnitude diagram for pho- 
tometrically selected GC/UCD candidates (grey dots, from 
iFaifer et alj|201lh , and spectroscopically confirmed GCs/UCDs 
(orange circles, this work ) . The three UCDs discussed in 
iNorris fc Kannappanl fcOllf ) are indicated as the orange stars. 
Several other objects are likely UCDs, but lie outside the region 
covered by HST imaging. This makes it impossible to definitively 
prove that their half light radii meets the necessary limit to be 
classified as UCDs (R c > 6-8pc). Upper Panel: Histogram of 
GC/UCD g' -i' Q colours for all GC/UCD candidates (grey his- 
togram, left axis), and spectroscopically confirmed GCs/UCDs 
(orange histogram, right axis). 
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Cote et al.l l200ll; iForbes et al.l |2004| ; iTamura et ail 
lHairidl2009l ; iForte et al.ll201ll ). The red GC distribu- 



tion is also similar to that of the diffuse light of NGC 3923, 
again a common fi nding, at least for the innermost re- 
gions of early-types jForbes et alj|2004 ITamura et al]|2006l ; 
iLee et all I200I iBassino et alj|200ct ). The slight remaining 
offset in slope between the red GCs and the diffuse light 
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Figure 4. GC surface density profiles. Solid dark lines in each 
panel are the galaxy V band surface brightness profile. Solid 
coloured lines in each panel (grey, blue and red) show the GC 
surface density profiles for the full GC sample, the blue subpop- 
ulation and the red subpopulation respectively. The profiles have 
been normalised to match the galaxy starlight profile at 60 arcscc- 
onds. The number counts of GCs within this region are expected 
to s uffer significant in completeness due to the high background 
Csee lFaifer et al.l201ll) . Shaded regions show the 1-cr poisson error 
regions. The purple dot-dashed line in each plot shows the best 
fit to the GC surface density profiles for GCs with galactocentric 
radii > 60 arcseconds. 



of NGC 3923 can most likely be explained as being due to 
interloping blue GCs that act to reduce the red GC slope. 

The effect of contamination of one GC subpopulation 
by the other has n ot been discussed in great detail in the 
literature. However, iForte et al.l (|2007l ) do discuss the possi- 
bility of contamination of the blue population between the 
blue peak and colour valley, by red GCs, in the galaxies 
NGC 1399 and NGC 4486. 

The mixing of blue GCs into the red sample due to 
the blue tilt i s also likely at work in the observation by 
ISikkema et al.l (|2006h that the colour profile of the red GCs 
is similar to that of the diffuse light of NGC 3923, but is 
offset toward the blue by around 0.05 mag. This int erpreta- 
tion is supported by the results of lFaifer et al.l (|201ll ) , where 
the use of Gemini/GMOS g and i photometry allows a bet- 
ter separation of the two populations, resulting in a red GC 
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colour profile that exactly matches that of the NGC 3923 
galaxy halo light. Such a close correspondence between the 
colours of the red GC s and the underlyin g integrated light 
is not unexpected. In iNorris et al.l (|2008T ) the stellar pop- 
ulation parameters of the diffuse light of NGC 3923 (at 
large radii) are found to be indistinguishable from those 
of the red GCs. This result is consistent with the case of 
the edge-on SO NGC 3115, where the stellar populations of 
the bulge component at ~2R C (measured along the minor 
axis), are likewise indis tinguishable from those of the red 
GCs jNorris et alj|2006h . 

3.3 Kinematic Analysis 

Figure [5] shows a histogram of the velocities of all 79 con- 
firmed GCs/UCDs, and for the blue (34), red (45), and 
brightest (15 objects with My < -10.5) subpopulations sep- 
arately, divided as described in Section 13.11 As can be 
seen from Figure \Sj unlike in our studies of NGC 4649 
^Bridges et alj|200rj ) or NGC 3379 (|Pierce et alj|2006t l. the 
recessional velocity of NGC 3923 (~1800 kms -1 ) is suf- 
ficiently high that contamination from Milky Way stars 
can be ignored. An examination of the environment around 
NGC 3923 also demonstrates that contamination of the sam- 
ple by GCs associated with nearby neighbour galaxies is 
likely to be negligible. Therefore it is assumed that all 79 
objects with velocities between 1100 and 2500 kms -1 are 
GCs/UCDs associated with NGC 3923. 

We choose to include in ou r sample the three UCDs dis - 
covered around NGC 3923 by INorris fc Kannappar] (120111 ) 
(including NGC3923-UCD3 which we spectroscopically con- 
firm here). We do this because the preponderance of the 
evidence is that these UCDs formed alongside the GCs of 
NGC 3923, and are merely the most massive members of the 
NGC 3923 GC system. The evidence for this is compelling; 
the UCDs luminosities and frequencies are consistent with 
being members of the NGC 3923 GCLF, the multiplicity 
of UCDs (at least 3, with 2 further high probability candi- 
dates just within the coverage of our GMOS imaging) makes 
a stripping origin unlikely, and most importantly their stel- 
lar populations are indistinguishable from those of the GCs 
of NGC 3923 |Norris fc Kannappanl l201ll . Norris et al. in 
prep). 

Overplotted on Figure [5] are the mean velocity and ve- 
locity dispersion of the various GC samples. The mean ve- 
locities of the blue and red subsamples differ by 25 kms -1 , 
however, this difference is not statistically significant given 
the 54 kms -1 uncertainty, and likely reflects small number 
statistics. Likewise, while the blue GCs appear to have a 
higher dispersion than the reds (the difference is 72 ± 40 
kms -1 ), this result is also not significant given the present 
sample size. The observation that the two GC subpop- 
ulations have similar, or only marginally different veloc- 
ity dispersion, is comm on to many studie s of early-type 
GC kinematics (see e.g. Beaslev et al.l |2004| ; iRichtler et al.l 
l2004l ; IWoodlev et al .112010 )7 Nevertheless, if increased statis- 
tics were to show that the blue GCs do indeed display a 
higher velo city dispersion this wo uld not be u nexpected, 
NGC 1399 JSchuberth et al.ll2010l ). N GC 4472 dZepf et all 
l200d ; IC6te et al.ll2003l i and NGC 4636 (|Park et alj|2010l ) all 
have blue GC subpopulations which are kinematically hotter 
than the red subpopulations. A higher velocity dispersion for 
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Figure 5. Histogram of GC velocities, upper panel shows all 
79 GCs/UCDs studied here, next panel shows the 45 velocities 
of the red objects, the second from bottom panel is for the 35 
blue objects. The lowest panel shows the velocities of all objects 
which have My < -10.5, roughly the magnitude limit of the o nset 
of UCD-like behaviour (see e.g. INorris fc~ Karmappan 201jj). Of 
these 15 objects, 3 are confirmed UCDs, 2 are high probability 
UCD candidates (being partially resolved in our GMOS imaging) 
and several others are likely candidates requiring HST size deter- 
minations for confirmation. The vertical da shed line in all pan els 
denotes the systemic velocity (1832 kms -1 [Carter et al,lll998f) of 
NGC 3923. 



the blue clusters could be explained naturally if some of the 
blue population are actually intra-cluster/intra- group GCs, 
or were accreted from dwarf galaxies stripped by the current 
host. 

Most interesting is the behaviour of the most luminous 
NGC 3923 GCs and UCDs. The bottom panel of Figure 
[5] demonstrates that unlike the case for the full sample, the 
most luminous objects show no preference for velocities close 
to the systemic velocity (si milar to the behavio ur seen in the 
M87 GC/UCD system bv lStrader et alj|201ll ). This obser- 
vation can have two possible explanations: 1) That the lumi- 
nous GC population is significantly contaminated by objects 
which are not NGC 3923 GCs, such as UCDs formed by the 
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stripping of dwarf galaxy nuclei. 2) That the most luminous 
GCs have orbits which differ from less luminous GCs, per- 
haps due to the preferential destruction through dynamical 
friction of massive (and more extended) GCs on orbits which 
bring them close to the centre of NGC 3923. Given the small 
numbers currently on hand it is hard to definitively prove 
which explanation is the correct one. However, we believe 
that 1) is unlikely, given the fact that it would likely require 
a significant fraction of the luminous NGC 3923 GCs to be 
interlopers. 

In the following analysis, NGC 3923 is assumed to 
ha ve an effective r a dius ( R c ) of 53.4 arcsec, as reported 
by Fukazawa et al.l (|2006T ) and originally determined by 
iFaber et all (jl989l ) in the B band. The use of this value, 
which differs from the value o f 43.8 arcsec o nds ( as mea- 
sured in the J-band) quoted in iNorris et al l (|2008l ). allows 
consistency between t he present study and the X-ray one of 
iFukazawa et all (|200rj) . Therefore, using this value of R c our 
last measured tracer is located at ~6.9 R c . The measurement 
of GC kinematics at such large radii is important, because 
at these radii the mass of the stellar component should be 
negligible. This makes the detection of a DM component 
much simpler than in the inner regions of galaxies, where 
the baryonic component dominates. 

Figure [6] displays histograms that highlight the distri- 
butions of azimuthal angle and galactocentric radius of the 
79 GCs. As expected from examination of Figure [3] there 
is a slight dearth of kinematically confirmed GCs located 
directly along the major axis of NGC 3923 (located at PA 
= 48°). It appears that this deficit is due to the combina- 
tion of a slight under density of candida te GCs in this region 
(as seen in Fig. 1 of iFaifer et al.ll201ll ) combined with the 
previously discussed problems of placing slitlets due to over- 
laps. The group of objects located at around 0° would pro- 
duce spectra which overlapped with any placed at around 
50°. The histogram of galactocentric distance demonstrates 
that velocities have been measured for GCs located at large 
galactocentric radii, up to 370 arcsec. 

3.4 Rotation of the GC system of NGC 3923 

The rotation of the GC system of NG C 3923 has been ex - 
amined using the procedure outlined in lNorris et al.l (|2008l ). 
A non-linear least squares fit was carried out to the equation: 



V(6) = V n 



5>o)+Vo 



(1) 



where V ro t is the amplitude of the GC rotation, is the az- 
imuthal angle, 8q is the position angle of the line of nodes, 
and Vo is the recessional velocity of NGC 3923. This proce- 
dure was carried out for each of the three GC samples; total, 
red subpopulation, and blue subpopulation. In practice the 
value of Vo was allowed to vary between samples. However 
the overall effect of this was small, as the mean velocities 
of the three groups varied by less than 10 kms -1 . The GC 
mean velocity (1802 ±17 km s -1 ) and systemic v elocity of 
NGC 3923 as determined by ICarter et~aH (|l998T > (1832 ± 
3 kms -1 ) differ by 30 kms -1 . This means that using the 
GC mean velocity instead of the systemic velocity could be 
masking some rotation. To assess whether or not this effect is 
significant, the rotation of the total sample was re-estimated 
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Figure 6. Upper Panel : Histogram of azimuthal an- 
gle (measured North-through-East and folded about 180°). 
Lower Panel : Histogram of galactocentric radius for the 79 
NGC 3923 GCs. The major axis position angle of NGC 3923 is 
48°. 



using the mean GC velocity ±30 kms - , in neither case was 
the result significantly changed. Figure [7] shows the best fit 
rotation curves derived using this procedure, in all three 
cases no significant detection of rotation is made. 

The lack of rotation observed in the GC system is con- 
sistent with the similarly small amount of rotation found for 
the integrated light of NGC 3923, either found using veloci- 
tie s determined from t he deep MOS observations presented 
in INorris et al.l (l2008h . or in smaller radii longslit studies 
by Carter et al.l (1 19981) (wh ich extended to 40 arcsec) and 
iKoprolin fc Zeilingerl feoOCt ) (which extended to 25 arcsec). 



A lack of significant rotation is also a commonly observed 
feature of the GC systems of large early-type galaxies, as 
discussed in Section [T] 

The lack of any significant differences in the observed 
kinematics of the GC subpopulations, and the high proba- 
bility of contamination of the "red" GCs by blue subpopula- 
tions members lead us to focus on the total GC population 
in all that follows. 



3.5 Velocity Dispersion 

Figure [5] presents the radial profile of the velocity dis- 
persion of the NGC 3923 diffuse light (green squares and 
stars) and total GC system (black line). The velocity dis- 
per sion of the GCs was estimated using the lowess estima- 
tor (|Gebhardt et al.lll994l ; lGebhardt fc Fischerlll995l ). which 
estimates the velocity variance at each datapoint by fit- 
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Figure 7. Velocity versus azimuthal angle for all 79 NGC 3923 
GCs. Blue circles are blue GCs, red squares are red GCs. The 
black dashed line represents the best-fitting rotation curve for the 
full sample, the blue and red dashed lines are the best-fitting rota- 
tion curves for the red and blue subpopulations. The implied ro- 
tation amplitudes are: 32.6±44.0 krns" 1 (red), 105.7±73.3 kms" 1 
(blue) and 40.4±37.9 kms -1 (total), none are statistically signifi- 
cant. The major axis position angle of NGC 3923 is 48°. 



ting a least squares linear fit to the individual points in 
a kernel (essentially in radius) around the point in ques- 
tion. Each datapoint in the fit is weighted by the inverse 
square of the distance between the point being fitted and 
the others within the kernel. The velocity dispersion is sim- 
ply estimated as the square root of the variance measured 
by the lowess method. As a consistency check, this method 
was compared with the velocity dispersion estimations pro- 
duced by a maximum likelihood method; for any reasonable 
bin size (minimum number of GCs per bin > 15), the max- 
imum likelihood method recovered velocity dispersions did 
not differ significantly from those found using the lowess 
method. 

The observation that the velocity dispersion profile of 
the NGC 3923 GCs is approximately constant, or slightly 
rising at larger radii, is a general finding in studies of the 
GCs of early type galaxies, as was discussed at length in 
Section [T] 

Of particular interest in Figure [H] is the fact that the 
velocity dispersions of the integrated light of NGC 3923 de- 
termined from MP S slitlets using the method outlined in 
iNorris et al.l (|2008l ) (green stars) are consistent with the ve- 
locity dispersions of the GC system at the same radii. The 
implications of these results are best understood through 
the examination of simple dynamical models of NGC 3923. 



3.6 Dynamical Models 

We use the axisymmetric, or b it-bas e d mo de ling algorithm 
descri b ed in iGebhardt et all (|200Cl. |2003l) . IThomas et ail 
l|2004l . |2005h and ISiopis et all |2009r i. These models are 



based on the method of ISchwarzschild (1 19791). an d 
lar model s are p resented in iRichstone fe Tremainej 
Rix et all (ll997l) . ICretton et all (|l999l ). IValluri et al] 



Cappellari et all (|2006l ) 



1984), 



2004). 



The modelling procedure consists of the following: (i) 
the surface-brightness distribution is deprojected into a 3d 
luminosity profile; (ii) a trial gravitational potential is cal- 
culated including the contribution from the stars and a dark 
matter halo; (iii) a large set of orbits is run in the trial po- 
tential; (iv) an orbit superposition model of the galaxy in the 
trial potential is constructed that matches the kinematical 
data subject to the constraints given by the 3d light pro- 
file and the trial potential; (v) steps (ii)-(iv) are repeated 
for a large variety of trial mass profiles. A x 2 - anai ysis then 
determines the best-fit model and its uncertainties. 

Due to the difference in the slopes of the profiles (sur- 
face brightness for stars and number for GCs), we model 
the stars and GCs separately, and then sum their \ 2 as a 
measure of the quality of the fit for the given mass profile. 
Thus, the only addition to the above procedure is that, in 
step (iv), we fit a different set of orbital weights for stars 
and clusters. Each model is also constrained by the surface 
brightness/number density profile of the particular tracer. 
In this case we use the observed V-band surface brightness 
profile for the diffuse light of NGC 3923, and the observed 
total GC sample number density profile as shown in Figure|4] 
and discussed in Section [4] Both stars and GCs are assumed 
to share the same flattening of e = 0.4, consistent with the 
E4-5 clas sification and the mea sured ellipticity of NGC 3923 
found bv lSikkema et all (|2007l ). 

The NGC 3923 models use a polar geometric layout as 
follows. We use 20 radial bins and 5 angular bins for the 
spatial sampling, and 15 velocity bins. The radial spatial 
binning ranges from 1 arcsec in width at small radii up to 
50 arcsec in the outer bins. The 5 angular spatial bins are 
equally spaced in sin(i), where i is measured from the major 
axis to the pole (i.e., angular bins are smaller closer to the 
major axis). The velocity bins are 100 kms' 1 wide. The 
average number of orbits per model is 40,000. The orbit al 
sampling follows the design in lThomas et al. I (|2004l2005l ). 

We fit for three parameters that determine the mass 
profile; these are the stellar M/L assumed constant with 
radius, the dark halo scale radius, and the dark halo density 
at R = 0. To represent the dark matter profile, we use a 
spherical cored logarithmic profile given as: 



pDM(r) 



3r, 



2 +r 2 



4ttG (r? +r 2 ) 2 



(2) 



where v c is circular ve locity at the core radius, r c 

|Binnev fc Tremaind[l987t) . 

IThomas et al] j2005l ). lForestell fc GebhardtJ (|2010l ) , and 
iMurphv et al] ( 201ll ) find that a cored logarithmic profile 
provides a better fit or at least as good a fit as an NFW 
profile. Given that we do not have kinematic data at large 
radii in large enough quantity to explore the various dark 
halo parameterizations, we only consider the one model. In 
all models, we include a black hole of mass 1 x 10 s Mq , which 
is the mass derived f rom its correlation with velocity disper- 
sion (|Gultekinll2Q10l ); the value of the black hole mass will 
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Figure 8. Upper Panel: Velocity dispersion versus radius for NGC 3923 GCs. Blue circles and red squares are the individual blue and 
red GC velocities, the black solid line shows the GC dispersion profile estimated using the lowess method described in Section I3.5I the 
orange shaded region shows the l-a uncertainties on the GC velocity dispersion profile. The blue and white line shows the dynamical 
model fit to the GC velocity dispersion profile, as described in Section 13.61 Lower Panel: This panel provides the comparable diffuse 
light measureme nts to those provid ed in the upper panel. The green squares are the velocity dispersions measured for the inner regions 
of NGC 3923 bv lCarter et al.l il998T), the gr e en sta rs represents the velocity dispersion measurements of the spheroid of NGC 3923 made 
using the method presented in lNorris et al ] {20081 '). The solid green line shows the result of fitting a dynamical model to the NGC 3923 
diffuse light as described in Section 13.61 The velocity dispersion profile of the GCs is also replotted for comparison purposes. As can be 
seen, the velocity dispersion of the diffuse light and GCs are identical over the common range of around 100 to 200 arcsec. 



not affect the dark matter profile and we only include it for 
completeness. 

The results of this fitting procedure are displayed in 
Figures [8l l9l and [TOl The best fit M/Ly in the central region 
is 8, id enti cal to that found for NG C 4649 bv lBridges et al.l 
(|2006h and lShen fc GebhardtJ teOlOT), and dose to the M/L v 
of 6 found for NGC 3379 by IPierce et"aH l|2006h . Mass-to- 
light ratios as high as these, while large, are only moderately 
inconsistent with the predictions of stellar population syn- 
thesis models for such old stellar populations (as the star 
symbols in the right panel of Figure [lOl demonstrate) . 

In the upper panel of Figure[5]the white and blue dashed 
line displays the best-fit model GC velocity dispersion pro- 
file. The fact that a model constrained by both the diffuse 
light velocity dispersion, and the GC kinematics (plus their 
respective number density profiles), is still able to success- 
fully reproduce the observed GC velocity dispersion profile 
is reassuring. 

In addition to determining the best-fit mass profile, 
orbit-based models provide a measure of the internal or- 
bital structure. Figure [9] presents the ratio of the radial to 
tangential internal velocity dispersion. We present this ratio 
for both the diffuse light (solid black line), and GCs (solid 
red line) of NGC 3923. In general, the large uncertainties 
prevent any strong conclusions on the orbital behaviour of 



the GCs and diffuse light. However, to first order the pro- 
files of both appear to be mildly radially biased at all but 
the largest radii. Additionally, the GCs and diffuse light of 
NGC 3923 appear to show similar behaviour over the region 
of overlap. Beyond a radius of around 130" both the GCs 
and diffuse light appear to become more isotropic (cr I /ot=l)- 

In the left panel of Figure [T0]the inferred enclosed mass 
profiles are displayed. Of immediate interest is the fact that 
the mass profil e dete rmined from X-ray observations by 
iFukazawa et al.l (|2006l ) is around a factor of two less than 
our total dynamical mass profile, though with similar slope 
at larger radii. In fact, the X-ray mass profile is lower than 
even our stellar only model (blue line). Therefore the dis- 
agreement cannot be due to differences in the derived DM 
halos of the models. However, the stellar mass and the X-ray 
mass profile could be made to agree, if a less extreme M/Lv 
of around 5 was chosen (which would also be more consis- 
tent with stellar population models), though this would then 
be inconsistent with the observed velocity dispersion of the 
integrated light and GCs of NGC 3923. 

Discrepancies of similar magnitude b etween X-rays 
and GCs have b een found in NGC 4636 (|johnson et al.l 
l2009h. NGC 1407 (iRomanowskv et al.ll2009r ), and NGC 4649 
( Shen fc Gebhardtll2010i ) . In fact disagreements between dif- 
ferent dynamical tracers within the same galaxy are a gen- 
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eral result of studies of galaxy dynamics. In particular, PNe 
commonly indicate the presence of little or no DM in galax- 
ies in which oth er tracers detect significant quantities e.g. in 
NGC 821 PNe (|Romanowskv et alj|2003 ; iTeodorescu et all 
l201(fl indicate little or no DM present, while stellar kin emat- 
ics (|Weiimans et al.|[2009l ; iForestell fc Gebhardtll2010l ) indi- 



cate between 13 and 18% of the mass of NGC 821 is DM 
within R e . As all dynamical tracers have their own strengths 
and weaknesses, only by combining several different tracers 
within a galaxy can we hope to achieve robust results. 

The best-fit parameters to our cored logarithmic profile 
are v c =450 ± 50 kms -1 , and r c =ll ± 4 kpc. The best-fit 
model has x 2 =22. For the stellar kinematics, we only use the 
dispersions measured in 13 circular apertures at various radii 
from 1.5"to 200". As discussed previously, we do not find 
significant rotation. For the globular clusters we use three 
spatial bins, and for each bin we consider both the velocity 
and velocity dispersion. Thus, the total degrees of freedom 
is 19. The reduced x 2 of 1.16 suggest that the dynamical 
model well represents the data. For the uncertainties, we 
rely on those models within a set A\ 2 value. We report 
68% confidence values, where we have marginalized over the 
other parameters; thus, the 1-sigma values in this case come 
from A\ 2 = 1. The confidence bands in Fig. 1101 represent 
the range in models that are within this \ 2 - 

Our best fit cored logarithmic model returns an enclosed 
dark matter fraction of 17.5^45% within R c , which rises to 
41.2j^gg% within 2R e . At the limit of our kinematic cover- 
age (370 arcsec, 6.9 R e ) dark matter comprises 75.6jljg' 8 % 
of the mass of NGC 3923. This distribution of dark mat- 
ter is similar to that measured in other massive early types 
using stellar absorption line kinemati cs or dynamical trac - 



ers, e.g. M87 (17.2% DM within R B iMurphv et al.ll201~lT) 
NGC 821 (13 - 18% DM w ithin R„. IWeiimans et al.ll2009l ; 
Forestell fc Gebhardtl201Ch . NGC 4594 (19% DM within R e , 
Bridges et all 120071 ). or 17 early- type galaxies in the Coma 



cluster (10-50% DM within R R . iThomas et all 12007] '). It is 
however slightly inconsistent with other studies which have 
found larger DM fractions using either S AURON observa- 
tions (median DM fraction 30% within Re. lCappellari et al.l 
I2006T I or o bservations of stron g gravitation lenses (38% DM 
within R„. |BoIton et alj|2008h . 

The right panel of Figure [TO] shows the behaviour of 
the M/L v profile of NGC 3923. At the limit of the binned 
data (269 arcsec) the inferred M/Lv of 25.9 is inconsis- 
tent with the constant M/Ly (of 8) at the 3.3cr level. This 
again confirms the existence of significant amounts of DM 
in NGC 3923. 



4 DISCUSSION 

An examination of the GC radial velocities has not de- 
tected any significant rotation in the NGC 3923 GC system. 
This finding means that the GCs display similar behaviour 
to the integrated l ight of NGC 3923 which also has little 



or no rotation (e.g 
20001: iNorris et all 



Carter et al. Hl998l ; lKoprolin fc Zeilingerl 
20081) . This lack of rotation in the GC 



systems of massive early-type galaxies is not uncommon 
with little or no rotation obse rved in the GC systems of 
NGC 4649 dBridges et al.ll2006l ). NGC 1399 dRichtler et al ' 



|2004| ; ISchuberth et alj l2010t ). NGC 3379 jBergond et al 
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Figure 9. The ratio of the radial to tangential velocity dispersion 
as determined using our dynamical modelling procedure outlined 
in Section I3.6I The solid black line shows the behaviour of the 
diffuse stellar light of NGC 3923, with the orange shaded region 
displaying the ±lcr errors. The solid red line shows the behaviour 
of the GCs of NGC 3923, with the dashed red lines denoting the 
±1<t errors. 



l2006f ). NGC 4472 dZepf et alj|2000l : ICote etHl |2003| ) . and 
NGC 4594 (|Bridges et al.ll2007l ). Some GC systems how- 
ever do display signifi cant rotation, for examp l e NG C 5128 
jWoodlev et all 120071) and M87 JCote et al.l l200ll ). The 
cause of this variation may be partly due to projection ef- 
fects. Nevertheless the current observations do not fit com- 
fortably with theoretical simulations of galaxy formation. 
These simulations tend to predict that significant amounts 
of angular momentum should be found in the outer regions 
of galaxies, especially in cases where the galaxies formed in 
major merger events in the relatively more recent past (i.e. 
z < 3. IVitvitska et alj|2002l ). 

To explain this disagreement it may be possible to call 
on angular momentum transport to move the angular mo- 
mentum beyond the radii to which current GC studies are 
presently sensitive; only an extension of these studies to even 
larger galactocentric radii would be able to determine if this 
is indeed the case. This would most likely prove impossible 
for all but the most massive galaxies with the richest GC 
systems, as smaller galaxies simply would not have enough 
GCs at sufficiently large radii to make such a measurement 
possible. 

Turning to the velocity dispersions of the GCs and in- 
tegrated light of NGC 3923 reveals that these also display 
indistinguishable behaviours over the range where they can 
both be measured (see Fig. [SJ. It is also intriguing that 
the value of the observed constant velocity dispersion of 
~200 kms -1 is indistinguishable from the di spersion mea- 
sured for the NGC 3923 group as a whole by iBrough et al.l 
(|2006h . using the redshifts of all 30 group members. It should 
be noted however, that there are indications that the in- 
tegrated light velocity dispersion may in fact be closer to 
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Fig ure 10. Lef t Pan el: Enclosed mass profiles, for the dynamical mass (red line), stellar mass (blue line), and X-ray mass profile from 
Rika zawa et al, I i2006h . The vertical dashed line shows the limit of the binned GC velocity data. Right Panel: Mass-to-light ratio profile. 
The apparent turn up of the model within 10 arcsec is due to the inclusion in the model of a nominal black hole of mass 1x10 s Mq. 
The star symbols re present SSP model predictions for M/Lv for the stellar populatio n given the age and metallicity measured for the 
population eithe r in lNorris et al.l l|2008h (larger radii point) or bv lThomas et al ] l|2005l l (small radii point). The SSP models come from 
iMaras ton (2005) for Salpeter IMF (upper/blue star symbols) and Kroupa IMF (lower/green star symbols). 



that of the red GC subsample, rather than the system as a 
whole. The weighted mean of the integrated light velocity 
dispersion for the four largest radii points is 198 ± 9 kms -1 , 
consistent with the 200 ± 22 kms" 1 found for the red GCs 
(Figure[5|, but marginally inconsistent with the velocity dis- 
persion of the blue GCs of 272 ± 34 kms . Thi s is compat- 
ible with th e halo-GCs connection discussed in lForte et al.l 
(|200ll2009l ). in the sense that most of the galaxy halo light 
is contributed by the diffuse stellar population associated 
with these red subpopulation clusters. Nonetheless this is 
not to say that the kinematics of red GCs alone would be a 
preferable tracer of the mass profile of NGC 3923. In prin- 
ciple, with a sufficient sample size, and hence the correct 
number profile and orbital information red and blue GCs 
should give the same mass profile. 

Because of the difficulty of measuring stellar absorp- 
tion line kinematics to large radii in early type galaxies, 
very few studies have been able to directly compare GC 
and i ntegrated light k i nema tics over the same radii. How- 
ever. ISchuberth et all (|2010l ) found that the red GCs and 
integrated light of NGC 1399 likewise display similar ve- 
locity dispersions over their common range of study. If this 
finding proves to be a general result of studies of GC kine- 
matics it will provide a very important clue to the for- 
mation history of galaxies. Such a close correspondence 
in kinematics, combined with the observed similarities in 
distribution (see Figure [4] and Section 13.21) . colours (see 
iForbes fc Fortd^OOllJsikkema et al.ll20od ). and stellar pop- 
ulations (see lNorris et al.ll2008l ). indicate a common forma- 
tion for the red GCs and diffuse light of NGC 3923. 

The observation that the velocity dispersions of the GCs 
and integrated light of NGC 3923 are constant at larger 
radii, despite the rapidly declining number density/surface 
brightness, is already strong evidence for the existence of 
significant amounts of DM in this galaxy. Our orbit-based 



models confirm that NGC 3923 has dark matter fractions 
consistent with those found in other studies of early type 
galaxies using integrated light or GC dynamical tracers. 

The observation that our estimates of dynamical mass 
(and hence M/Ly) are la rger than those made w ith X-ray 
observations is intriguing. IShen fc Gebhardtl (|2010f ) find dis- 
crepancies of similar magnitude and in the same sens e in 
their study of NGC 4649, as do iMurphv et ail (|201ll ) for 
their study of the st ellar kinematics of M87. In contrast, 
I Johnson et"aH (|2009l ) find disagreement in both senses in 
the same galaxy, NGC 4636, where within 10 kpc they find 
the GCs indicate ~twice as much mass as the X-rays, be- 
tween 10 and 30 kpc GCs and X-rays agree, and then beyond 
30 k pc X-rays indicate 4-5 tim es as much mass as the GCs. 
Like IShen fc Gebhardtl (|2010f ) we suggest that the observed 
discrepancy may indicate problems with the X-ray mod- 
elling. Perhaps other forms of non-thermal pressure support, 
such as magnetic fields, microturbulence, or cosmic-rays are 
present in the X-ray gas, but are not yet adequately treated 
in the X-ray mass determination. Additionally, the assump- 
tion that the X-ray emitting gas is in hydrostatic equilibrium 
may need to be relaxed. This is a particular concern in dis- 
turbed galaxies such as NGC 3923, where its obvious shell 
structures may imply a relatively recent merger event. These 
results act to reinforce the fact that only by combining as 
many mass tracers as possible within individual galaxies is 
it possible to overcome each tracers shortcomings and arrive 
at a robust estimate of galaxy dynamics and DM content. 



5 CONCLUSIONS 

This paper has presented Gemini/GMOS spectra of 79 GCs 
associated with NGC 3923 measured using MOS and N&S 
techniques. Their radial velocities have been combined with 
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GC number densities profiles, stellar integrated light kine- 
matics, and stellar surface brightness profiles and used to 
constrain the DM content of the NGC 3923. Our main con- 
clusions are: 

(i) There is no significant evidence for rotation in any of 
the GC sub-samples (total, blue, or red) out to the largest 
projected radius measured (6.9R C ). 

(ii ) In common with other studies of early type galaxies 
(e.g. ISchuberth et all l20ld ; iPark et al.1 120101 ) the blue GC 
subsample appears to have a higher dispersion (272 ± 34 
kms" 1 ), than the red GC subsample (200 ± 22 kms" 1 ). 
But with our current sample size this difference is not yet 
statistically significant. 

(iii) The GC velocity dispersion profile for the total GC 
sample is consistent with being constant with radius. 

(iv) The velocity dispersion of the integrated light of 
NGC 3923 at large radii is also constant, and is consistent 
with the GC velocity dispersion profile at the same radii. 

(v) We find some evidence that the diffuse light and GCs 
of NGC 3923 have radially biased orbits within ~130". 

(vi) By applying axisymmetric orbit-based models to 
the GC and integrated light data we find that NGC 3923 
is 17.5±I;|% DM within R c , A1.2t^%% within 2R C , and 
75.61 within 6.9R C . These DM mass fractions for 
a group elliptical are indistinguishable from those found 
using similar methods for ellipticals in environments from 
isolated (NGC82I) to cD (M87), perhaps indicating the 
universal nature of the DM profile. 

(vii) The total dynamical mass within the radius 
of our last dynamical tracer (6.9R C ) is found to be 
1.5t^ 5 xl0 12 Mq. 

(viii) The inferred GC+integrated light dynamical mass 
and X-ray mass profiles are discrepant. The X-ray mass pro- 
file is around 50% of that of the dynamical mass. 
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Table 1. Confirmed NGC 3923 GCs. ID, RA, Dec, X, Y, g, r, i, velocity and velo city errors for all 79 confirmed NGC 3923 GCs/UCDs. 
The ID refers to the GCs position in the photometry table of lFaifer et al.l j201llV The X and Y positions refer to po sitions relative to 
the galaxy centre. Name denotes the ID given to objects studied in iNorrisetal lf [|2008h or[N orris fc Kannappar] j201 if ). 



ID 



RA 



Dec. 



X 





(J2000) 


(J2000) 


(arcmin 


26 


177. 


70500 


-28 


78370 


2.74 


35 


177. 


70690 


-28 


82370 


2.64 


39 


177. 


70700 


-28. 


81640 


2.63 


55 


177. 


78160 


-28. 


.77420 


-1.29 


58 


177. 


70860 


-28 


83960 


2.55 


90 


177. 


71090 


-28 


78090 


2.43 


103 


177. 


77170 


-28. 


81650 


-0.77 


128 


177. 


78560 


-28 


77520 


-1.50 


158 


177. 


78340 


-28 


81080 


-1.38 


162 


177. 


78300 


-28. 


76930 


-1.36 


176 


177. 


78150 


-28. 


.77910 


-1.28 


193 


177. 


77960 


-28 


83760 


-1.18 


206 


177. 


77880 


-28 


82970 


-1.14 


214 


177. 


74260 


-28 


84840 


0.76 


243 


177. 


73280 


-28. 


80520 


1.28 


251 


177. 


77530 


-28. 


80570 


-0.96 


275 


177. 


77380 


-28 


83190 


-0.88 


277 


177. 


77380 


-28 


82810 


-0.88 


312 


177. 


77010 


-28 


76440 


-0.69 


318 


177. 


76990 


-28. 


78640 


-0.67 


336 


177. 


76700 


-28. 


80560 


-0.52 


344 


177. 


76680 


-28 


76360 


-0.51 


365 


177. 


76470 


-28. 


83240 


-0.40 


375 


177. 


76410 


-28. 


79420 


-0.37 


380 


177. 


77950 


-28. 


82690 


-1.18 


396 


177. 


76120 


-28 


78020 


-0.22 


419 


177. 


75850 


-28 


82620 


-0.07 


453 


177. 


75500 


-28. 


76550 


0.11 


458 


177. 


75470 


-28. 


.81330 


0.13 


484 


177. 


75260 


-28 


76210 


0.24 


192 


177. 


75160 


-28 


78010 


0.29 


502 


177. 


75060 


-28. 


80160 


0.34 


503 


177. 


75060 


-28 


78290 


0.34 


530 
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72510 


-28 


80540 


1.68 


534 


177. 


74800 


-28. 


80790 


0.48 


535 


177. 


74780 


-28. 


83600 


0.49 


552 


177. 


74640 


-28. 


.76710 


0.56 


557 


177. 


74590 


-28 


83100 


0.59 


593 


177. 


74050 


-28 


79270 


0.87 


595 


177. 


74140 


-28. 


77960 


0.82 


601 


177. 


74040 


-28. 


77930 


0.88 


617 


177. 


73830 


-28 


80010 


0.99 


621 


177. 


73840 


-28 


83810 


0.98 


ooy 


177. 
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-28. 


76740 




683 


177. 


73260 


-28 


77300 


1.29 
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177. 


73240 


-28 


78010 


1.30 


713 


177. 


73050 


-28 


77580 


1.40 


725 


177. 


72900 


-28. 


.80110 


1.48 


760 


177. 


69530 


-28. 


82620 


3.25 


779 


177. 


72420 


-28 


79050 


1.73 


791 


177. 


72330 


-28 


76960 


1.78 


824 


177. 


79280 


-28. 


80270 


-1.88 


847 


177. 


68510 


-28. 


85610 


3.78 


859 


177. 


71780 


-28 


83340 


2.06 


872 


177. 


71640 


-28. 


80990 


2.14 


903 


177. 


71360 


-28. 


84480 


2.29 


940 


177. 


80100 


-28. 


80420 


-2.31 


952 


177. 


78790 


-28. 


80630 


-1.62 



Y 

(arcmin) 



So 
(mag) 



'0 

(mag) 



'0 

(mag) 



V 

(kms ) 



Name 



22 


.52 


22.05 


21.85 


1586.5 


± 


66.0 


NSB228 


21 


.43 


20.81 


20.53 


1590.8 


± 


35.9 




22 


.35 


21.68 


21.38 


1726.6 


± 


28.9 


NSB197 


21 


.31 


20.78 


20.59 


2295.2 


± 


47.3 


NSB65 


22 


.68 


22.15 


21.98 


2203.2 


± 


74.2 


NSB253 


23 


.08 


22.61 


22.35 


1851.2 


± 


69.2 


NSB394 


20 


61 


20.04 


19.80 


2307.9 


± 


34.8 


NGC3923-UCD3 


21 


.76 


21.27 


21.07 


2144.1 


± 


45.0 


NSB104 


22 


.69 


21.99 


21.63 


1809.0 


± 


31.7 




22 


12 


21.52 


21.23 


1810.9 


± 


26.9 


NSB167 


21 


22 


20.60 


20.32 


1995.9 


± 


50.4 




21 


80 


21.13 


20.78 


2046.0 


± 


32.5 


NSB108 


22 


.73 


22.20 


21.99 


1752.7 


± 


45.4 




22 


.86 


22.13 


21.81 


1925.2 


± 


55.0 


NSB333 


19 


.99 


19.42 


19.17 


1478.0 


± 


28.6 


NGC3923-UCD2 


21 


.36 


20.60 


20.22 


1790.4 


± 


28.4 




22 


.88 


22.13 


21.69 


1638.9 


± 


40.9 


NSB332 


22 


.68 


22.02 


21.67 


1945.0 


± 


64.2 




21 


.81 


21.22 


20.88 


1555.1 


± 


35.7 


NSB110 


23 


.03 


22.39 


22.06 


1555.8 


± 


39.0 




19 


.53 


18.88 


18.58 


2114.7 


± 


29.9 


NGC3923-UCD1 


22 


19 


22.01 


21.82 


1805.4 


± 


53.1 


NSB221 


22 


52 


21.81 


21.47 


1936.0 


± 


56.0 


NSB225 


21 


.23 


20.52 


20.15 


1438.4 


± 


27.7 




21 


11 


20.97 


20.78 


1159.5 


± 


50.2 




22 


70 


22.10 


21.81 


1190.8 


± 


44.5 




21 


.88 


21.40 


21.20 


1564.0 


± 


60.2 




23 


16 


22.86 


22.42 


1618.9 


± 


61.2 


NSB513 


21 


.60 


21.12 


20.93 


1617.4 


± 


62.7 




23 


.26 


22.65 


22.37 


2114.4 


± 


38.9 


NSB450 


22 


.22 


21.55 


21.22 


1760.6 


± 


15.2 




21 


.56 


20.76 


20.36 


1627.5 


± 


26.6 




22 


.56 


21.96 


21.57 


1807.9 


± 


29.6 


NSB232 


21 


.28 


20.68 


20.42 


1764.7 


± 


32.0 


NSB64 


21 


.56 


20.88 


20.54 


2148.3 


± 


26.1 




22 


97 


22.27 


21.88 


1457.0 


± 


70.4 


NSB360 


23 


.33 


22.60 


22.26 


1749.1 


± 


74.1 


NSB492 


23 


.03 


22.35 


22.07 


1621.8 


± 


73.0 




21 


.87 


21.27 


20.99 


2114.6 


± 


52.7 




22 


,83 


22.10 


21.74 


1796.2 


± 


55.6 




22 


.79 


22.20 


21.93 


1666.4 


± 


42.5 


NSB311 


21 


.87 


21.18 


20.86 


1949.9 


± 


24.2 




22 


.61 


22.06 


21.86 


1840.7 


± 


48.2 


NSB247 


21 


.62 


21.03 


20.72 


1784.4 


± 


51.1 


NSB86 


22 


.80 


22.08 


21.67 


1755.7 


± 


53.9 




23 


.36 


22.75 


22.42 


1622.1 


± 


48.6 


NSB498 


22 


.83 


22.30 


22.08 


2038.4 


± 


88.5 


NSB322 


21 


72 


21.08 


20.77 


1664.7 


± 


24.9 


NSB99 


20 


.50 


19.89 


19.62 


1714.4 


± 


26.3 




21 


.42 


20.85 


20.61 


1861.3 


± 


43.8 




21 


.73 


21.14 


20.85 


1733.2 


± 


26.3 


NSB106 


21 


.62 


21.03 


20.79 


1756.4 


± 


42.5 




21 


.89 


21.45 


21.06 


2024.5 


± 


68.1 




21 


.66 


21.00 


20.68 


1970.9 


± 


59.7 




21 


.87 


21.27 


21.02 


1762.7 


± 


26.2 


NSB112 


22 


70 


21.98 


21.61 


1796.6 


± 


57.9 


NSB279 


21 


26 


20.74 


20.55 


1494.4 


± 


56.9 




21 


.66 


21.00 


20.67 


1738.2 


± 


25.3 


NSB93 



1.35 
-1.05 
-0.61 
1.92 
-2.00 
1.52 
-0.62 
1.86 
-0.27 
2.22 
1.63 
-1.88 
-1.41 
-2.53 
0.06 
0.03 
-1.54 
-1.31 
2.51 
1.19 
0.04 
2.56 
-1.57 
0.72 
-1.24 
1.56 
-1.20 
2.45 
-0.42 
2.65 
1.57 
0.28 
1.40 
0.05 
-0.10 
-1.78 
2.35 
-1.48 
0.81 
1.60 
1.62 
0.37 
-1.91 
2.33 
2.00 
1.57 
1.83 
0.31 
-1.20 
0.95 
2.20 
0.21 
-2.99 
-1.63 
-0.22 
-2.31 
0.12 
-0.00 
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Globular Clusters and the Dark Matter Content of NGC 3923 

Table 2. Tabic Q] continued. 



ID RA Dec. X Y g r' i' V Name 

(J2000) (J2000) (arcmin) (arcmin) (mag) (mag) (mag) (kms -1 ) 
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22 


,25 


21 


.50 


21 
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1988 
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± 


31 


4 


1356 


177.69350 


-28 
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3.34 


-0.47 


22 


,02 


21 


.36 


21 


.06 
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39 


4 
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-28 
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3.93 
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21 
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20 
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20 
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± 


23 


1 
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-28 
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4.53 
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22 
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21 
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